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The simple interpretation of PAO's UHECRs' penetration depth measurements suggests a transi- 
tion at the energy range 1.1 — 35 • 10 18 eV from protons to heavier nuclei. A detailed comparison of 
this data with air shower simulations reveals strong restrictions on the amount of light nuclei (pro- 
tons and He) in the observed flux. We find a robust upper bound on the observed proton fraction 
of the UHECRs flux and we rule out a composition dominated by protons and He. Acceleration 
and propagation effects lead to an observed composition that is different from the one at the source. 
Using a simple toy model that take into account these effects, we show that the observations requires 
an extreme metallicity at the sources with metals to protons mass ratio of 1:1, a ratio that is larger 
by a factor of a hundred than the solar abundance. This composition imposes an almost impossible 
constraint on all current astrophysical models for UHECRs accelerators. This may provide a first 
hint towards new physics that emerges at ~ 100 TeV and leads to a larger proton cross section at 
these energies. 



Among the most interesting results of the Pierre Auger 
Observatory (PAO) are the penetration depth measure- 
ments [lH3j: the observed mean depth where maxi- 
mal number of secondaries are generated, (X max )(E), 
and the fluctuations of this quantity RM S{X max ){E) = 
&(X max )). When compared with extensive shower sim- 
ulations [J-tZl both (X max ) and <j(X max ) are consistent 
with a protonic composition around 1 EeV . However, 
at higher energies, (X rnax ) and o(X max ) fall below the 
protonic simulated values, towards values estimates for 
intermediate mass (e.g. iron) nuclcQ . At first sight this 
suggests a transition, around 10 EeV, from protons to 
intermediate mass nuclei [lip. We explore here the im- 
plications of the penetration depth data combined with 
the observed energy spectrum lfj- Hf . 

Consider a cosmic ray flux composed of N species each 
with a fraction /,-, a mean penetration depth, (X max }j, 
and RMS variation Uj . At each energy these are related 
to the measured mean and RMS values as: 



N 



) fj (-^-max^j — (-^maxji 



(1) 



N 



fj (°j + (Xmax)j) ~ (X ma x) 2 ~ <j{X max ) 2 . (2) 

In the following we examine several possible solutions of 
these equations for different compositions. 



1 The interactions between UHECRs and atmospheric particles 
takes place at CM energies of ~ 100 TeV. This is about a hun- 
dred times higher than energies in which cross sections have been 
measured. 

2 While PAO has the best statistics, the High resolution fly's eye 
(HiRcs) and the telescope array (TA) observatories compo- 
sition results are consistent with proton dominated composition 
all the way up to the high end of the spectrum. This disagree- 
ment among the different observatories is still to be settled. 



Consider, first, a mixture of protons (denoted p) and 
another arbitrary component, denoted 0. This arbitrary 
component may be a single species or a combination of 
a few. The condition (Tg > yields an upper bound on 
f P (E): 
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Fig. Q] depicts the upper limit obtained using the ob- 
served values of (X max ) and <j(X max ) and the simulated 
values of (X max } p and a p . The maximal proton fraction 
is smaller than 50% at E > 10 19 eV and it decreases be- 
low 30% at higher energies. This upper bound depends 
only on the observed PAO data and on the shower simu- 
lation results for protons (it does not even depend on the 
shower simulations for nuclei). It is independent of the 
acceleration or propagation of the UHECRs. When we 
replace the hypothetical ingredient with any composition 
of He, N, Si and Fe the upper bound on f p at the two 
highest energy bins is smaller by « 10% than the upper 
limit derived using <j$ = 0. 

Protons and Fe survives best the interactions with the 
cosmic photon field while propagating from the sources 
to Earth [3] ■ As such they are the most natural UHECR 
ingredients. However, with just two components the 
system of eqs. CD)-© is overdetermincd. Wilk and 
Wlodarczyk [14[ have shown that there is no consistent 
solution for protons and Fe within the error bars of the 
PAO data and the shower simulations. 

Eqs. CD)-© nave a marginal (x 2 /dof ~ 1) solution 
for a mixture of protons and He (see also [15J). How- 
ever, when propagation effects are taken into account 
this composition can be ruled out. He nuclei photodisin- 
tegrate rapidly on their way from the sources to Earth. 
For simplicity we neglect redshift effects and assuming 
a uniform distribution of sources. Under this approxi- 
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FIG. 1. An upper bound on the proton fraction of the 
observed flux as a function of energy, calculated using the 
PAO data and different extensive air shower simulations: 
QGSJET01 H, QGSJET-II [| , Sibyll2.1 and EPOSvl.99 
[f| (see legend). 
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FIG. 2. The total observed total flux (blue dotted line with 
diamonds) and the calculated secondary proton flux arising 
from He disintegration (eq. J3}) (red dashed line with circles), 
for a composition with a maximal number of primary protons 
(satisfying eq. Q) and He (color online). At high energies 
the number of secondary protons exceeds the observed flux. 



mations the fraction of He nuclei with energy E surviv- 
ing and reaching earth is Fczk h<1 (E) = Ire /hi 1, 
where Ihe is the mean free path for He photodisin- 
tegration and lh is the horizon size. For a given ob- 
served He flux, (1 — f p (E))J(E), the flux at the source 
is (1 - f p (E))J(E)/F GZI<Hc (E). Since F GZKhf {E) « 1 
most of the He nuclei disintegrate producing secondary 
protons with energy E/A. The resulting secondary pro- 
ton flux, J PtSec (E) i|| 
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3 (i) The protons GZK distance is comparable to the horizon dis- 
tance at these energies, (ii) We have overestimated here the GZK 
distances of 3 He and 2 H as equal to the GZK distance of 4 He. 



Using the upper bound on f p {E) obtained earlier (eq. 
([3])) we find (see fig. [2]) that this secondary proton flux 
is larger than the maximal proton flux allowed. 

Kampert and Unger (l5l | have shown that a mixture 
of protons, He and intermediate elements like N, Si and 
Fe can provide a solution of eqs. flU and ^ for the ob- 
served composition which is within the uncertainties of 
the observed data and the simulations. However, as we 
have seen for p and He, to determine the sources' compo- 
sition we need to take propagation effects into account. 
Moreover, since different species are accelerated differ- 
ently within a given accelerator, acceleration should also 
be considered. To examine these effects we consider a 
toy model based on only two components: protons and 
Fe. As mentioned earlier, this composition cannot satisfy 
equations ([T]) and and other intermediate elements in 
addition to Fe are needed. Therefore, when we exam- 
ine a proton and Fe composition we consider only the 
overall spectra and (X max ) and we ignore, for simplic- 
ity, the RMS data. This simple example is sufficient for 
demonstrating the nature of the problem. 

Any electromagnetic acceleration process that acceler- 
ates protons to energy E accelerates nuclei (with charge 
Z) to energy □ ZE. This suggests a natural 

nana 

ex- 
planation for the transition in composition: the source 
accelerates protons to a power law energy distribution, 
E~ a up to some maximal energy where a gradual cutoff 
begins. The same source accelerates nuclei to the same 
power law but up to an energy that is Z times larger. 
This naturally produces a heavier observed composition 
than the one at the source and may suggest that the drop 
at very high energies in the UHECR flux is not necessarily 
due to a GZK effect but simply due to lack of available 
accelerators that can accelerate UHECRs to extremely 
large energies [18j |. We characterized the accelerator's 
cutoff by an unknown function of the rigidity, g(E/Z), 
with < g(E) < 1 and g(E) = 1 at low energies. 

However, before adopting this model propagation ef- 
fects should also be taken into account. Like acceleration, 
propagation in the IGM magnetic field depends on the 
rigidity. On the other hand GZK attenuation depends 
on the nucleus at hand. In this energy range (1.1-35 
EeV), it is negligible for protons but it is significant for 
all nuclei. We characterize the propagation effects us- 
ing Fqzk Fc = Ipe/lh [IB]- Under these assumptions the 
total observed UHECRs flux is: 



J(E) = c p g(E)E- a + c Fe F GZKpe (E)g{—) 
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4 In the diffusive shock acceleration, a proton and a nucleus, with 
an atomic weight A, crossing the shock front will have the same 
Lorentz factor, and therefore a nucleus will be A times more 
energetic. This suggest that we need to compare nuclei at energy 
E with protons with energy E/A and not E/Z. This will add 
a factor of (A/Z) a ~ l ~ 3 to the composition ratio. This factor 
doesn't change qualitatively the results 
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Where c p , cpe are normalization factors for the proton 
and Fe nuclei fluxes respectively. c p is obtained by the 
condition g(E) = 1 at the minimal energy. The normal- 
ization of cp e is such that cp e /c p equals the Fe nuclei to 
protons number ratio at the source (before acceleration 
Using the proton fraction f p (E) = c p g(E)E~ a / 'J(E) and 
g(E/26) = 1, which is valid over the relevant energy 
range (1.1-35 EeV) eq. ((5j becomes: 
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We solve eq. (pj for f p {E) using the measured {X max ) 
and the simulated values for protons and Fe. Now that all 
the quantities at the l.h.s of ([6]) are known we fit a power 
law (E~ a ) to the l.h.s to to obtain a and cp e . The best 
fit results are: a = 2.1 ±0.1 , c Fe /c p = (2.2 ± 0.6) • 10~ 2 , 
with x 2 /dof = 0.21. 
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FIG. 3. The l.h.s of eq. (J75J) ( with errorbars, blue circles 
for Fe, green squares for Si) and the best fit for the r.h.s of 
cq. © (green solid line for Si, blue dashed line for Fe, color 
online). The Si results are shifted to the right for clarity. The 
best ht results for Fe and protons: a = 2.1 ± 0.1 , cp e /c p — 
(2.2 ±0.6) ■ 10~ 2 , with xV do / = 0.21 and for Si and protons: 
a = 1.9 ± 0.1, csi/c p = (5 ± 1) • 10~ 2 , with y^jdoj = 0.4. 



The number ratio we obtained corresponds to a mass 
ratio of ~ 1 : 1. A similar analysis with Si instead of 
Fe yields a = 1.9 ± 0.1, c Si /c p {E) = (5 ± 1) • 10~ 2 , with 
X 2 /dof = 0.4. The mass ratio is again « 1 : 1. As 
both Fe and Si show almost the same trend a composite 
solution that will satisfy both the {X max ) and cr(X max ) 
data [TBI ] will have similar features and a ~ 1 : 1 mass 
ratio between the protons and the metals. 

Note that we have neglected secondary protons aris- 
ing from photodisintcgration. Taking those into account 
would have resulted in even higher metallicity. Shifting 



5 Within this acceleration model the number of Fe nuclei accel- 
erated to energies larger than 26E equals to cp e /c p times the 
number of protons accelerated to energies larger than E (Note 
that |22H use somewhat different notations) 



the observed energies by about 20%, as suggested by com- 
parison of the PAO spectra with the spectra observed in 
other main UHECRs observatories [21|, does not change 
qualitatively our results. Interestingly, this extremely 
heavy composition at the source is comparable to the up- 
per limit obtained using the angular distribution Q of 
these UHECRs. These ratios of N Fe (> 26E)/N p {> E) 
are > 0.072 for the VCV catalogue and > 0.084 for cor- 
relations with Cen A. 

Finally, we note that the spectral index, found here, is 
much harder than what observed in lower energies, a = 
3 [19(. This arises from the GZK attenuation affecting 
nuclei at these energies. Such hard spectra were obtained 
in detailed propagation simulations [li, 17, 2(J. 

We have shown that the PAO penetration depth mea- 
surements and the penetration depth numerical simula- 
tions yield a robust upper limit on the observed proton 
fraction of the UHECRs flux. This limit drops below 50% 
at energies higher than 10 EeV and below 30% at higher 
energies. These measurements are inconsistent with the 
composition of 75% protons and 24% He, that is common 
in the Universe and they require a significant fraction of 
intermediate mass nuclei. 

The conversion of the observed composition to the 
composition at the source depends on the acceleration 
and propagation. Using a simple toy model we find that 
the protons to metals mass ratio at the source should be 
about 1:1. This metallicity is larger by a factor of a 
hundred than the solar metal abundance of ~ 1%, which 
reflects typical metallicity in the Universe. This high 
metallically puts a new severe constraint on the sources, 
since objects dominated by nuclei heavier than He are 
rare in the astrophysical landscape. 

Active galactic nuclei (AGNs) are natural UHECRs 
accelerators [23], [24| • Most AGNs can accelerate metals 
to ~ 10 20 eV and protons to an order of magnitude lower. 
This would produce, naturally, the observed composition 



transition 



WL 17 1 as well as the cutoff in the spectrum 
at higher energies [HI]. However this would require a 
very heavy composition, whereas AGNs typically show 
solar- like mctallicities [25| . 

Gamma ray bursts (GRBs) are another natural UHE- 
CRs accelerators [H (see however (13, If UHECRs 
are produced by GRBs' internal shocks, they are com- 
posed of the original material of the jets. One has to in- 
voke, in this case, a very efficient nucleosynthesis within 
the jets and survival of the produced nuclei during the 
acceleration [2^, [3(J. It is not clear that this can be 
achieved generically. Recall that the UHECR output of 
GRBs should be comparable or larger than their 7-rays 
output. Alternatively if UHECRs are produced by ex- 
ternal shocks they will be composed by the circum-burst 
winds surrounding the star. GRBs' progenitors, Wolf - 
Rayet stars, have He, C and O dominated wind whose 
composition is too light. 

A variation on this theme was proposed by Liu and 
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Wang [3l| who invoke, instead of regular GRBs, low- 
luminosity GRBs which they describe as "semi relativis- 
tic core collapse supernovae that involve an engine activ- 
ity" and call "hypernovae" . These are also based on Wolf 
- Rayet progenitors. However, low-luminosity GRBs jets 
do not penetrate the stellar envelope and the observed 
emission is produced by a shock breakout [32|. Thus, 
it is not clear how could these bursts accelerate UHE- 
CRs in the first place. Furthermore, as sources of heavy 
UHECRs they suffer from all problems mentioned earlier 
concerning regular GRBs. 

A rapidly rotating young pulsar with a strong but rea- 
sonable magnetic field (~ 10 13 G) can accelerate Fc to 
UHECR energies (H. Fang et al. [H suggested that 
young pulsars are UHECR sources and the origin of the 
heavy composition is the Fe rich crust. However, X-ray 
observations of neutron stars suggest the existence of an 
atmosphere composed of proton and light elements above 
the crust (see [33( and citations therein). Acceleration of 
this atmospheric component will result in a light compo- 
sition. This poses a serious doubt concerning this model. 

Overall, while astrophysical heavy UHECRs sources 
cannot be ruled out with absolute certainty, the strict 
constraint on the composition obtained here (that con- 
strain both the protonic and the He components) makes 
the UHECRs sources puzzle even harder to solve. In 
particular it rules out the most natural sources, AGNs. 
There is no single clear model that naturally produces 
UHECRs with such a composition. One may search 
for acceleration processes that are not rigidity dependent 
and favor heavy nuclei over light ones, however such a 
process is not readily available. 

Given this situation one can consider the following al- 
ternatives to the heavy composition. First, the obser- 
vational data might be incorrect or somehow dominated 
by poor statistics: these results are based on about 1500 
events at the lowest energy bin and on only about 50 at 
the highest one. The possibility of a miscalculation in 
the shower simulations is unlikely. Different simulations 
0-01 obtain comparable results. However, the simula- 
tions depend on the extrapolations of the proton's cross 
section from the measured energies of a few TeV to the 
range of an UHECR - atmospheric nuclei collision, which 
are factor of 100 higher in energy. Is it possible that this 
extrapolation breaks down? A larger cross section than 
the one extrapolated can explain the shorter penetration 
depth. If so these findings might provide hints of a new 
physics that set in at energies of several dozen TeV [35 1. 
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